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Abstract Short laser pulse in wide range of wave- 
lengths, from infrared to X-ray, disturbs electron- ion equi- 
librium and rises pressure in a heated layer. The case 
where pulse duration tl is shorter than acoustic relax- 
ation time ts is considered in the paper. It is shown that 
this short pulse may cause thermomechanical phenom- 
ena such as spallative ablation regardless to wavelength. 
While the physics of electron-ion relaxation strongly de- 
pends on wavelength and various electron spectra of sub- 
stances: there are spectra with an energy gap in semi- 
conductors and dielectrics opposed to gapless continuous 
spectra in metals. The paper describes entire sequence 
of thermomechanical processes from expansion, nucle- 
ation, foaming, and nanostructuring to spallation with 
particular attention to spallation by X-ray pulse. 
PACS: 79.20.Ds, 65.40.De, 81.16.-c 



1 Introduction 

There are many industrial applications using short pulse 
lasers. New exciting possibilities are connected with de- 
velopment of the X-ray lasers. Ablations by the long and 
short pulses differ qualitatively. The first - evaporate, 
boil, and, at higher fluences, move matter by ablative 
pressure created in hot plasma corona. While the re- 
lease of pressurized layer is the main process in the case 
of short pulse. It is shown below that this is true for any 
laser wavelength. 

Let's consider short pulse. Irradiation with sufficient 
intensity transfers substance into warm dense matter 
state. In condensed phase, the cohesive properties are 
important. The cohesion is result of interatomic attrac- 
tion. Due to stiff behavior of solids and liquids, their 
expansion in rarefaction wave is very different from ex- 
pansion of gas. Stiff means that the typical moderate 
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expansions zAp/p ~ (0.1 — 0.2) cause the order of magni- 
tude pressure drops and even change of sign of pressure - 
from compressed to stretched state, here Ap is the den- 
sity drop due to expansion. Stretched metastable layer 
under tensile stress appears as a result of the release of 
the warm dense matter with its stiff response. 

In metastable state substances are sensitive to tem- 
perature and degree of stretching. Nucleation probabil- 
ity exponentially depends on amplitude of negative pres- 
sure. Therefore, sharp nucleation threshold appears on 
this amplitude. The threshold depends on temperature. 
Nucleation is followed by development of two-phase layer 
composed of condensed phase and voids. If the laser 
heated layer melts, then nucleation takes place in liq- 
uid. In this case, expansion of the two-phase layer leads 
to foaming. The foaming may cause formation of the 
surface nanorelief [1]. 

The metastability and nucleation are well-known for 
release initiated by shock coming to the surface from 
the bulk of a target. But in the case of shock sent by 
a long laser pulse, an ion beam, or by an explosion of 
chemical explosive, near nucleation threshold the two- 
phase layer locates far from the target surface. In the 
case of large scale, the development of the two-phase 
layer cannot disturb the surface of a target because the 
nucleation layer and the target surface are far from each 
other and are independent from each other. On the con- 
trary, the short pulse lasers initiate foaming very close 
to the surface - since the attenuation depth datt for 
X-ray photons, or thickness of a skin layer Sgkin^ may 
be as small as ten nanometers. In this case the foaming 
strongly interacts with the surface, eventually produc- 
ing frozen surface structures. Therefore we can use such 
terms as nanofoam, nanostructures, or nanospallation to 
describe the situation with small depth datt- 

Photon absorption and collisional processes are de- 
fined by photon energy and electronic structure. Infrared 
(IR) and visible radiation excite valent electrons whereas 
X-rays are absorbed mainly by internal shells in the one- 
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photon interactions. This is why the X-ray absorption is 
quahtatively similar for metals, semiconductors and di- 
electrics. Action of quanta from the one electron- Volt 
range of energy depends strongly on existence of the for- 
bidden gap A. In metals, they are absorbed mainly via 
inverse Bremsstrahlung in the skin layer. In cases with 
gap, the seed Keldysh ionization, inverse Bremsstrahlung 
heating of ionized electrons and electron avalanche con- 
trol the rise of number of conduction electrons Ueix^t) 
[SlElllllS] . Significant heating takes place when plasma 
frequency for these electrons oOpiirie) overcomes laser 
frequency uol during the laser pulse. During the rest of 
the pulse, the substance with a gap absorbs laser energy 
in a skin layer. This means that during the rest of the 
pulse the absorption becomes similar to the absorption 
of metals. This greatly increases spatial density J • cm ~^ 
of absorbed energy. Estimates show that, in these condi- 
tions, laser electric field strength is comparable to atomic 
fields of external electrons, their wave functions are dis- 
torted by electromagnetic wave, and probabilities of the 
multiphoton and tunnel ionizations are significant. 

Relaxation of electrons to equilibrium, after the end 
of a pulse, depends on a band structure. In metals, free 
electrons cool due to electron- ion energy transfer and due 
to thermal conductivity. The same is true for semicon- 
ductors which pass to metallic state during their melting. 
In substances which keep the gap after a pulse, the con- 
centration of free electrons Ue and their temperature Tg 
decrease after the end of a heating laser pulse as a result 
of recombination, diffusion of electrons and holes, and 
electron heat conduction. Three-body recombination is 
usually more significant than radiative recombination. 
Relaxation time teq is of the order of 1-10 ps for all 
these cases. 

An acoustic response time ts is necessary to decrease 
pressure few times in the laser-heated layer with thick- 
ness (It- It equals to ts = / Cs in bulk targets, or 
ts = df /cs in foils with small thickness df < dr^ where 
Cs is sound velocity. 

Electrons are light, their velocities are high and the 
electronic thermalization time Tg is small - typically 
it is at the femtosecond range. In metals, this corre- 
sponds to rather high values of Tg ~ 1 eV. This case 
is considered here. At room temperatures Tg ~ 300 K 
the e-e relaxation in metals is slow, since e-e collision 
frequency z^gg is small. If Tg ~ 1 — 10 fs then at the 
picosecond time scale we have two thermodynamic sub- 
systems: electrons and ions. Total pressure p = Pe-\-Pi is 
composed of partial pressures. The contribution Pe > 
because Tg — > 0; in the one-temperature state we 
have Pe = 0. 

In our conditions, typical absorbed energy is Eabs ^ 
(0.1 — l)Ecoh'', where Ecoh is heat of sublimation; e.g., 
for Al Ecoh ^ 3eV/atom. Then, after e-i relaxation, ion 
temperatures are in the kiloKelvin (kK) range: 

Ti ^ 4 {Eabs/ lev • atom-^) kK, 



Pi ^ 19(K6./leV-atom-^)(n/6xl0^^cm-^)GPa, (1) 

if heat capacity is ^ Sks, and Gruneisen parameter is 
Fi ~ 2. The parameter F = V{dp/dE)v links pressure 
rise and fast absorption of energy at the isochoric stage. 

At the two-temperature stage t < teq, electron tem- 
peratures Tg are much higher than T^. For metals, in 
the Fermi-gas approximation, we have Tg ~ y2E^^/7, 

Tg = 17Z-^/^(n/6xl0^^cm-^)^/^(T;g/leV-atom-^)^/\K 

for te = Fe/FF < 1, ksFF = Ef. Here Ee - Eabs 
is electron energy per atom, not per electron. Z is the 
number of electrons per ion, 7 = 7r^ng/c^/2T^F is the 
electron heat capacity constant written in the Fermi-gas 
approximation. 

Electronic contribution pe is significant when Tg ^ 
Fi, then Fe — Fi ^ Fe. In this case, in metals. 



Pe 



Pau ^ i^au J^au \ \ ^ ^ 



(2) 



where te is normalized temperature Tg, and Pau = 
2.94 X lO^GPa, Uau = 6.757 x 10^^ cm -^ 1/nJu = 
0.53A(Bohr radius), Eau = 27.2 eV are atomic units 
(au) for pressure, density, etc, 

EF/Eau = {l/2){3^^ Zn/Uauf^^ 

is Fermi energy. The interpolation ([2]) has the right limits 
for te < 1, tg > 1. 

Pe ^ 2.8 Z^/^(n/6 X 10^^cm-^)^/^(Tg/leV)^ GPa (3) 

if te < 1. 
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Fig. 1 (Color on line) Comparison of the ablation thresh- 
olds (for incident fluence) as function of Al and tl : three 
squares |[3j, filled circle ^23j, two empty circles - this work. 
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We call a laser pulse short if tl < tg. Usually the 
condition teq < ts holds. In such a case, the rarefaction 
release, driven by envelopes smaller mass than the 
Pi -release, and nucleation takes place in one-temperature 
state. In the case when teq > ts expansion, subsequent 
stretching and nucleation are connected with p^. In the 
two-temperature state, the vapor- solid and vapor-liquid 
coexistence curves are shifted in the direction toward the 
two-phase region. This is a result of a blow out by the 
electronic pressure of atomic system, composed of atoms 
which attract each other. At the same absorbed energy, 
the pressure Pe is few times smaller than p^, since the 
degenerate and classic gases of electrons are softer and 
their Gruneisen parameters = 2/3 is smaller than 
Fi ^ 1.5 — 2.5 for solids above the Debye temperature. 

The main condition for nanostructuring and spalla- 
tive ablation is tl < ts- Other parameters, such as 
wavelength (from IR to X-ray) and spectra {A = or 
Z\ 7^ 0), are less significant. Therefore, spallative abla- 
tion is an important mechanism for removal of material. 
For IR and visible quanta hv ^ 1 eV acting on metals 
and semiconductors, theory of spallative ablation is well 
developed [6l[7l[8]. The theory is developed to the point 
where fine and secondary consequences of the metastable 
decay such as nanostructuring in case of a wide focal spot 
dspot dT are predicted (HE]. 

For hv ^ leV and transparent dielectrics, the New- 
ton rings [TOt fTT] were never observed in spite of sev- 
eral attempts p!2l[T3l fT4] . It seems that interplay of two 
circumstances makes this observation difficult. Indeed, 
significant absorption Eabs together with small absorp- 
tion depth, similar to skin depth in metals, are achieved 
above threshold fiuence Fi)rkd for optical breakdown. 
At the same time, heating Eabs is a very sharp func- 
tion of absorbed fiuence Fabs near this threshold [21 
[SlIlllS]. Therefore, spallative ablation is limited to the 
narrow region near Fbrkd-, because the value of heat- 
ing is restricted Eabs < Eum-up, Eum-up ^ ^•'^Ecoh 
[H]. Above this limit, cohesive property is weak against 
strong stretching in a hydrodynamic rarefaction wave - 
expansion proceeds similar to expansion of heated gas 
(that is, without spallative plate). The cohesive prop- 
erty is responsible for creation of a spallative plate and 
spallative cupola p!6]. The cupola is necessary for the 
interference which results in appearance of the varying 
in time Newton rings [16]. 

The second circumstance is connected with the width 
of the gap A. The semiconductors such as Si and GaAs, 
with rather narrow Z\, metallize during melting. In the 
meanwhile molten dielectrics remain in dielectric state. 
In this case, the cupola is dielectric, and Newton inter- 
ference oscillations are weak (weak oscillations due to 
presence of oxide film have been detected in [13j). 

The Newton rings are bright manifestation of exis- 
tence of spallative ablation. Appearance of rings means 
that a light wave interferes between spallative cupola and 
the rest of the target. This is a very surprising example 



of spallative plate so thin that it is even transparent (!) 
to light - the skin depth is of the order of 10 nm. In 
more customary spallation by a long laser pulse [iTt fTS] , 
the plate is much thicker. 

Contrary to the case with quanta hv ~ 1 eV and 
non-transparent substances like metals and semiconduc- 
tors, the thermo mechanical effects for X-ray irradiation 
are not investigated. This is the first attempt to study 
thermomechanics, metastability and nucleation induced 
by X-rays. This subject is interesting in connection with 
fast progress in developing of X-ray lasers and in con- 
nection with a whole number of experimental papers on 

The paper is organized as follows: First, we show how 
the ratio r^/ts infiuences the maximum pressure cre- 
ated by absorption of laser energy. Then, new results 
concerning freezing of nanostructures at a late stage are 
presented. After that, the theoretical model of spalla- 
tive X-ray ablation, and experimental findings connected 
with this model, are described. It is shown that as a re- 
sult of the conditions, presented in the following lines 
(i), (ii) and (iii), the threshold Fau for the X-ray spalla- 
tive ablation is extremely low, in comparison with other 
cases with different laser wavelength \l and durations 
Tl- This is illustrated in Fig. [TJ The three mentioned 
conditions are: (i) 100% absorption of X-rays (no re- 
ffection); (ii) negligible diffusion and heat conduction 
loses out from the small attenuation depth datu ^ind 
(iii) smaller energy densities necessary for spallative ab- 
lation, in comparison with evaporative ablation of equal 
amount of material. 



2 Duration of a laser pulse and amplitude of an 
acoustic response 

Fast laser pulse transfers matter into energy containing 
state similar to a state of a chemical explosive behind 
a front of a detonation wave. This transfer is a base for 
hydrodynamic release, metastable decay and spallation. 
As was said, the pulse is short if its duration tl is com- 
parable or shorter than acoustic time tg. In this case, 
a pathway at a thermodynamic phase plane consists of 
two parts. One is passed during a pulse r^, while the 
other - during an acoustic response ts- The first part 
corresponds to the heating along the isochor p = piniUai 
with complications concerning the two-temperature de- 
tails. The second part is formed by an approximately 
isentropic release along an adiabatic curve. It intersects 
the coexisting curve and penetrates into the two-phase 
region. 

The states of material near the target surface, irra- 
diated by a long pulse tl ^ ts, are located near the 
coexistence curve. The hotter state relates to the larger 
absorbed laser intensity. Pressure created by a long pulse 
equals approximately to the pressure Psat-vap of satu- 
rated vapor. The later is limited by pressure Pcr in the 
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critical point, e.g., for Al pcr ^ 0.4 GPa. The critical 
pressure is significantly below pressures ([T]), (|3]) achieved 
during spallative ablation. 

For the same final temperatures Tshort\bin and Tiong\bir 
of material at the coexisting curve, in cases with short 
and long pulses, the maximum pressure Pshort\max at 
the short pulse thermodynamic pathway is much higher; 
here the subscript "bin" marks the coexistence curve also 
called the binodal. A short pulse pathway deviates sig- 
nificantly from the coexistence curve into the high pres- 
sure condensed phase region. This is why the pressure 
Pshortlmax IS higher. The Tshort\bin corresponds to in- 
tersection of the part 2 of the short pulse pathway with 
the coexistence curve. The first paper, where the path- 
way has been used, was the paper pjj. Now this useful 
conception is widespread p!|[6l[8l[24] . 

The temperature Tshort\bin\abi corresponding to spalla- 
tive ablation threshold is significantly below the criti- 
cal temperature Tcr- Increase of the absorbed fiuence 
Fabs above Fau, increases Tshort\bin above Tbinlabi- 
There are distinct fiuence F^v and temperature Ti)in\ev 
above which the spallative layer disappears [H] . The ra- 
tio Tcr/T})in\ev depends on material properties and de- 
tails of the two-temperature stage. It seems that it is 
larger for Au in comparison with Al. 

For F > Fev^ material expands without spallation 
plate. This regime is different from the "phase explosion" 
[25] by long pulse with Tiong\bin ^ ^cr- Release of high 
pressure p > 0, created isochorically by short pulse, 
produces expansion with high rate of stretching du/dx^ 
where u and x are along expansion direction. After nu- 
cleation in metastable state inside two-phase region the 
inertia of expanding matter infiates bubbles. The inertia 
is significant since the rate du/dx is high. Whereas in 
the case of the "phase explosion", the expansion of two- 
phase mixture is driven more slowly by weaker forces 
connected with pressure difference ~ Psat-vap —Pout be- 
tween pressure inside bubbles and pressure outside the 
target surface. 

Importance of the ratio TL/tg is illustrated in Fig.O 
The acoustic time for this case is tg ~ 100 nm/ (5.4 km/s) = 
20 ps. We see that positive and negative pressures begin 
to decrease in their amplitude when the ratio tl / ts be- 
comes larger than unity. Distributions of pressure in the 
bulk Al target are shown at the instant t = 10 ps after 
the maximum of intensity / oc exp(— t^/r£) of pump 
pulse. We use two-temperature hydrodynamic code de- 
scribed in Ref. t26j. Relaxation time teq is of the order 
of 3 ps. 



3 Freezing of nanostructures 

As previously stated, a short pulse initiates a sequence 
of stages: (i) release, (ii) metastable state, (iii) nucle- 
ation, (iv) evolution of foam. Nucleation takes place if 
absorbed fiuence is above nucleation threshold F^uci • It 
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Fig. 2 (Color on line) Strong decrease of maximum pressure 
for long pulses. 



is important that nucleation does not mean tacitly that 
spallative plate will run away. The spallative ablation 
threshold Fau and Fnuci are separated. The separa- 
tion {Fabi - Fnuci)/Fabi bclougs to the few % range. 
As was shown in [Ij, during release, the foam infiates to 
thickness comparable to the distance between the foam 
and free surface of a target. Therefore it perturbs surface 
causing appearance of the surface nanorelief. 

Development of the nanorelief is a dynamical phe- 
nomenon, caused by deceleration of free surface, infia- 
tion of foam, and surface tension resistance to infiation. 
There is the appearance time tappr when it develops. 
This time is much larger than tg because fiuence F 
is ^ Fabu and the nanorelief appears near the stop- 
ping point of free surface. Expansion velocities at this 
stage are small, and the nanoreleif develops slowly. For- 
mation of the nanorelief may be experimentally observed 
as changes in refiectivity. There are changes due to sur- 
face nanorelief and due to absorption in thick foam under 
surface [27j. 

In Al, and in many other materials, the ablation 
threshold Fau is higher than the melting threshold Fmeit , 
and then at F > Fnuci bubbles appears inside the 
molten layer. There are two cases with slow development 
of foam: one when Fnuci < Fabs < Fau and another 
when Fabs is slightly above Fau- In the first case, the 
foam remains closed under the surface. An example of 
this is shown in Fig.O In the second case, there is process 
of slow detachment of the spallative plate. The plate is 
connected by the random net of the liquid filaments with 
foam. The filaments are stretched and break off one af- 
ter another during detachment of the plate. During this 
process, a nanobrush from a forest of standing filaments 
appears. 

The subsequent evolution is material and target (bulk 
target versus foil, spot radial size) dependent. Let's com- 
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ment first the dependence on material properties. In 
heavy metals late stages of foam development are espe- 
cially slow. This is significant because there is competi- 
tion between the late links of the foam development cas- 
cade and the rate of conductive cooling. E.g., density and 
material strength ratios for Au vs. Al are 19.3/2.7 and 
20/13 [7J, and gold foam moves slower. In Au, the layer 
between foam and melting front is thinner than in Al, 
freezing temperature is higher 1337/933, molten layer 
is thicker, and thermal conductivity 318/237 is higher; 
while heat capacity and surface temperature 2.5 kK at 
ablation threshold are approximately equal. 

Cooling due to electron heat conduction is important. 
In liquid state, the foam under target surface F < Fau 
and nanobrush F > Fau finally disappears because the 
surface tension collapses bubbles with low vapor pressure 
in the foam and smooths down the threads and surface 
bumps. Bubbles collapses when supporting them ten- 
sile stress decreases to zero. Contrary to this smoothing, 
the undersurface bubbles F < Fau and "nano-brush" 
F > Fabi remain in the final relief if conductive loses 
are fast enough to freeze them. The nano-brush may be 
accompanied with frozen bubbles under. Comparison of 
Au vs. Al shows that in the case of Au, this is easier. 
This picture describes appearance of black gold in re- 
cent experiments p!|[28] . Focusing of the laser beam and 
thickness of foil are also significant, since tight focusing 
improves cooling, while in thinner foils the reservoir to 
adopt heat from molten layer is smaller. 

Molecular dynamic (MD) simulation of development 
of foam and its freezing requires huge computer power 
and smart multi-processors algorithm, because foam oc- 
cupies large volume and inflation/freezing processes are 
very slow. In the work |lj the first part concerning the 
inflation of foam has been done. Here, we present new 
results describing freezing. 

In order to include the electron thermal conductivity 
in the MD code, the electrons as classical particles were 
added to atom subsystem. Each "electron" is assumed 
to associate with its host atom/ion. An electron always 
allocates the position of host, but it has its own veloc- 
ity. The neighboring atoms may exchange their electrons 
with some frequency, determined from known experi- 
mental thermal conductivity of liquid aluminum/gold at 
melting point. The pair of atoms, where electron-electron 
exchange takes place, is chosen randomly from the list 
of neighbors of each atom. The list contains neighboring 
atoms within cutoff distance. 

In addition to the frequency of electron-electron ex- 
change between atoms, we adjusted the mass of "elec- 
trons" , in order to obtain the characteristic time of electron- 
ion energy exchange in collisions. The electron- ion scat- 
tering depends on the velocity vectors of an electron with 
mass mi as a projectile and target ion with mass m2 
according to equations: 

vi = (mi vi + m2 V2 + m2vn) / (mi + m2) (4) 




Fig. 3 (Color on line) Freezing of Al foam as result of con- 
ductive cooling. Green and red colors correspond to solid and 
liquid Al. The melt ing/recrystallizat ion front moves very slow 
at the time interval to which the left picture t = 154 ps be- 
longs. At the right picture t = 468 ps the overcooled liquid 
layer with T = 570 K is partly frozen, the underlying bubbles 
are strongly deformed but survived. 

V2 {mivi + m2V2 - mivn)/{mi + m2) (5) 

where v is an electron- ion relative speed, and and 
V2 are the velocities of particles in the laboratory system 
after collision. The unknown unit vector n is assumed 
to be distributed isotropically and is produced by using a 
generator of random directions. The Eqs. (|4l5p conserve 
the momentum and total energy of electron and ion sub- 
systems. Moreover, such a combined Monte- Car lo-MD 
approach guarantees conservation of local charge neu- 
trality. 

Our tests indicate, that the described above electron- 
electron exchange and electron-ion collision procedures 
lead to energy diffusion through dynamic net of atoms, 
simulated by MD, and give a correct solution of contin- 
uum heat conduction equation for liquid aluminum/gold. 

Two stages of evolution of Al nanobubble chain are 
shown in Fig. [3l Fluence in this case is slightly below 
spallative ablation threshold Fau- Therefore spallative 
plate keeps its connection to target after the process of 
inflation of bubbles and freezing. Our approach differs 
from the approach applied in the important paper [29], 
where the one-dimensional (ID) finite-difference scheme 
has been used. The scheme takes into account ID ther- 
mal conductivity and is solved parallel to MD simu- 
lation. The ID scheme employs transversally averaged 
temperature fields from the parallel MD simulation. 

It is important that in case of foam our method allows 
to describe spatial separation of thermal fluxes pumping 
heat through threads and walls with small cross-section. 
This greatly delays cooling through foam. Therefore the 
base of the foam, which is located from the bulk side, 
cools faster than the film between the foam and free sur- 
face as it is shown in Fig. [3] (b). A liquid layer above 
the bubbles in Fig. [3] (a) has temperature gradient from 
1380 K on its free surface to 1340 K at its internal sur- 
face contacting with bubbles. The liquid layer shown in 
Fig. [3] (b) is in supercooled state, comp. with [30j. Its 
temperature is T 570 K which is below melting point 
933.6 K on 360 K, and temperature gradient across the 
layer is small in this case. Separated nanocrystallites 
start to grow from two free surfaces of the liquid layer. 
As [3] indicates a bubble position is not random but cor- 
related with positions of second and third neighbors at 
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least. The freezing process fixes this correlation with 
formation of surface nanostructures having correlation 
length a few inter-bubble distance. 

The foam may be partially frozen from the bulk side, 
when, in the case with fluence above ablation thresh- 
old Fabs > Fahu still liquid film between the foam and 
free surface of target begins to transform to the spalla- 
tive plate and the process of breaking of the liquid- 
solid threads starts. In this case the forest of the frozen 
threads appears (black gold [UEH]). MD simulation shows 
that the gold foam has a thick multi-level structure with 
significant drop of temperature and expansion rate du/dx 
through levels from one level to another toward the bulk 
side, see movies |31J. Detailed description of the MD sim- 
ulation results, concerning foaming and freezing, is be- 
yond the scope of our paper. Here we have to show that 
a short X-ray laser pulse may cause similar thermome- 
chanical phenomena as an optical pulse. Thick foam is 
formed if the fluence F is significantly above the thresh- 
old Fahi' The "flashes" of the successive nucleations dur- 
ing a process of formation of thick foam take place when 
positive pressure p is lowered by the decrement Pstr 
[27] . where Pstr is material strength. At acoustic stage 
pressure is halved from the initial value to the value cor- 
responding to the established compression wave 0. If 
Pstr does not depend on temperature then the distances 
between successive flashes are equal. In case PstriT) the 
distances is growing toward the bulk since temperature 
decreases and PstriT) increases. 

Near threshold Fau the foam evolution time ~ 0.1 — 
1 ns greatly exceeds collisional relaxation times ~ teq ^ 
1 — 10 ps for hv 1 eV and X-ray lasers. Therefore 
the picture with melting and foam may be applied to 
the case of dielectrics heated by short pulse X-ray laser. 
Similar phenomena may appear in case of metals and 
semiconductors. In the case of LiF considered below the 
ablation threshold Fau is lower than the melting thresh- 
old Fmeit' The work p!5] is devoted to comparison of 
the spallative ablations from liquid and solid states. If 
we rise fluence up to Fmeit in LiF then deceleration of 
free surface to almost stopping point will be not possi- 
ble, because to slow down motion it is necessary to fulfill 
the condition F ^ Fau- This means that in LiF near 
Fahi there are no melting/freezing processes described 
above. But other substances, e.g. Al, have Fmeit < Fau- 
In these cases X-ray laser will cause nanostructuring. 



4 Theoretical model of X-ray action 

Dependence of the attenuation depth dattihv) taken 
from [32] is shown in Fig.Hl The right arrow corresponds 
to the Ag X-ray laser described in next Section. The 
depth datt is large for high energy photons. Then the 
acoustic response ts = datt/cs corresponding to hard 
photons is long. For energy Kuol = 12 keV, which will 
be achieved soon at XFEL ^33j, it is tg ~ 50 ns. This is 
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Fig. 4 Variation of datt{hiy) with photon energy hu. 

very large time in comparison with the picosecond time 
scale of the atomic processes. 

The estimate of the spallative ablation threshold for 
these hard photons is 



Fabl = C datt riat Fc, 



200 J /era 



where ( = 0.2 — 0.4 is a coefficient [7lfT5] in relation 
Eabi = C^cohi Fahi is energy per atom at the thresh- 
old, Uat is the atom concentration in solid state. Future 
XFEL lasers [33] wih have O.lkeV < Tiujl < 12keV pho- 
tons, ultrashort durations tl ^ 20 fs, and fluence up to 
several hundred J/cm^. Irradiation by such pulse can 
cause spallative ablation of huge sub-millimeter piece of 
condensed target. 

Equations describing the electron- ion non-equilibrium 
stage and hydrodynamic motion are 

dx/dt = pdx = p^dx^ ^ p^du/dt = —dp/dx^^ (6) 



P 



dt 



dqe 



-p. 



^dEi/p _ dqi 



dt 



dx"" 



du 

dx^ p 
du 



Pi 



dx° 



(7) 



(8) 



(9) 



dt Ox"" Ui2 

This is system of equations in Lagrangian coordinate x^ 
corresponding to an initial position of a material point. 
The initial density profile inside the target is = 

— oo) = p^. Values X, P = Pe^Pi^ Ef^'^ = neUi2 + 
Ee^ Ei^ Tie in (|6]|9]) are functions on variables x^, t; Ue 
is electron concentration ng in the conduction band; 
£^ei = <^ (Tg — Ti) is an electron-ion energy exchange 
rate; ^e,i = —{p^e,i/p^)dTe^i/dx^ are heat fluxes; j = 

— {pD/p^)dne/dx^ is diffusion flux of electrons; t\jQ i are 
electron and ion thermal conductivities. 

Primary electrons and primary holes are produced 
during X-ray pulse acting on LiF. Their kinetic and 
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potential energies are Ee\prim ~ 10 — 30 eV and un, 
Ee\prim + Un = hv = 89.3 eV. These energies are mea- 
sured from the bottom of the conduction band. Primary 
electrons and holes are non-equilibrium ones - they are 
described by the distribution function - they can not 
be described by temperature Tg and concentration ng. 
The non-equilibrium electrons and holes relax to ther- 
malized state with secondary electrons and holes through 
Auger processes, impact ionization, and three-body re- 
combination. Relaxation time is Trei ~ 1 ps. The value 
Trei is smaller than the pulse duration tl = 7 ps. There- 
fore, if we exclude the short time lapse Trei in the be- 
ginning of the laser pulse, then we can include contribu- 
tion of primary particles into electron energy budget ([7j) 
through energy conservation (radiative loses are small). 
It is supposed in equations ([7j) and ([9]) that the laser 
source {p^ / p)Q supplies energy directly to secondary 
electrons and produces secondary electrons at the rate 
Q/ui2, where Ui2 ^ A is an ionization potential of sec- 
ondary electrons, Z\ ^ 14eV is width of forbidden gap 
in LiF. 

Every primary electron produces {hu = 89.3eV) / {ui2-\- 
Ee) ^ 6 secondary electrons. Solution of system (|6]l9]) for 
parameters = 7ps, F = lOmJ/cm^, LiF, datt = 
28 nm gives ne\max = (1 - '^%)nat, Te\max ^ 2eV, 
where \max corresponds to the maximum values. These 
values are achieved at the end of X-ray laser pulse. The 
parameters correspond to the experiment described in 
next Section. Electrons are classical since their kinetic 
energy (3/2)/c^Te ~ 3eV is larger than Fermi energy 
Ep ^ 0.5 eV corresponding to our case with low concen- 
tration of free electrons. 

The energy transfer rate is calculated as E^i = a (Tg — 
Ti) ^ aTe = AEe^ A = 2a/3kB^ since heating of 
a lattice by electrons is significant only when electrons 
are much hotter than a lattice. Lattice temperature for 
F = lOmJ/cm^ is Tat\max ^ TOOK. Pressure wave 
initiated by this fast increase of temperature is shown 
in Fig. [51 The wave travels to the bulk side of a target 
(to the right side in Figure) from the irradiated surface, 
while temperature profile remains "frozen" into matter. 
The point x = is an initial position of the surface. The 
profile of the wave contains tail with negative pressure 
Pneg (tcnsilc strcss). The amplitude \pneg{t)\ is enough 
to cause spallative ablation. Crater depth is 40-50 nm. 
These findings agree with experimental results described 
below. 



5 X-ray experiment 

The experiment has been performed with the Ne-like 
Ag soft x-ray laser (XRL) facility at JAEA Kansai Pho- 
ton Science Institute, working at transient collisional 
scheme [22 1 I34 1 I35] . The XRL beam with an energy ~ 
1 microJ and the horizontal and vertical divergences of 
12 mrad x 5 mrad, respectively, was focused on a LiF 




100 



Fig. 5 (Color on line) Pressure and temperature profiles in 
LiF irradiated by X-ray laser pulse. 

crystal of 2 mm thickness and 20 mm diameter, by using 
a spherical Mo/Si multilayer mirror of 1050 mm radius 
of curvature (see Fig. [6]). The total energy on the LiF 
crystal of the XRL beam after passing 200 nm Zr filter 
and reflecting from the focusing mirror was ~ 170 nanoJ 
in a single shot. The luminescence of stable color centers 
(CCs) [36, 37, 38^139] formed by XRL radiation, was used 
to measure the intensity distribution in the XRL laser 
focal spot [221 135]. After irradiation of the LiF crystal 
with the XRL, the photo-luminescence patterns from the 
color centers (CCs) in LiF (shown in Fig. [6]) were ob- 
served by using a confocal fluorescence laser microscope 
(OLYMPUS model FV300). An OLYMPUS BX60 mi- 
croscope in visible differential mode and an atomic force 
microscope (AFM, TOPOMETRIX Explorer), operated 
in the tapping mode, have been used for measurements 
of the size of ablative spot. As it was shown in our previ- 
ous experiments [35], only about 6% of full laser energy 
is concentrated to the best focus spot of ~ 200 square 
microns. This corresponds to energy ~ 5mJ/cm^ or 
laser intensity ~ 7 x 10^ W/cm^. 

Two types of experimental investigations of XRL ab- 
lation threshold of LiF crystals were done. In the first ex- 
periments the Zr filter has been removed, and the XRL 
beam expands without attenuation by the filter. This 
powerful beam has been focused on the surface of LiF 
crystal. In Fig. [Tithe AFM image of the focal spot of this 
XRL beam is presented. The image is obtained after a 
single laser shot. Ablation of crystal is clearly seen at the 
AFM image and traces. The ablation threshold for LiF 
irradiated by a single shot is 10.2mJ/cm^. From the 
traces in Fig. [71 we could see that the ablation depths 
varied between 30 and 55 nm. These values are close to 
the theoretical crater depths calculated above. 

In the second type of experiments, the Zr filter was 
settled inside the propagation path of the XRL beam. 
Three shots have been done in the same focusing spot 



8 



Inogamov, Zhakhovsky, Faenov et al. 




jt 1 t t 1 t t t 



-3 -2 -1 +1 +2 +3 +4 

Fig. 6 (Color online) (a) The experimental set up for recording the Ag XRL beam patterns near the best focus position on a 
LiF crystal, (b) Sketch of motion of a LiF crystal during experiments, (c) The patterns of XRL beam focusing spots recorded 
on a LiF crystal at -3 mm to +4 mm from the best focus position. 




Fig. 7 (Color on line) AFM image and traces, taken through orthogonal directions, of the ablative spot on LiF crystal, 
irradiated by single shot with full (10.2 mJ/cm2) laser intensity of XRL beam pulse. 



with fiuence of 5 mJ/cm ^ for each shot. We could see in 
Fig. [8] that in this case^ a crater, with an ablation depth 
of about 50 nm, appears on the surface of the crystal. 
This is similar to the first experiment with single more 
intensive XRL shot. 



The ablation threshold obtained in our experiments 
is much smaller in comparison with previous experiments, 
see Fig. [H Our threshold is 3400, 300, and 10 times 
smaller than thresholds for nanosecond and femtosec- 
ond Tiisapphire lasers, and for nanosecond 46.9 nm soft 
XRL, respectively. 



6 Conclusion 

It is shown that short pulse of XRL causes thermome- 
chanical response as in cases with optical lasers. There 
is spallative ablation as result of appearance of tensile 
stress which overcomes material strength above abla- 
tion threshold. Near threshold stretching of melt may be 
accompanied by nanostructuring. Nanostructures freeze 
down if conductive cooling is fast. Value of the X-ray 
threshold is small in comparison with irradiation by longer 
wavelengths and/or longer pulse. 

Work has been supported by the RFBR grant No. 
09-08-00969-a (NAI, VVZh, VAK, and YuVP). This re- 
search has been partially supported by the Japan Min- 
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1 laser shots E = 5 rnJ/crn^, S.SxlO^ W/cm^ 3 laser shots E = 5 mJ/cm^, S-SxlO^ W/cm^ 




Visible image Visible image AFM image 



Fig. 8 (Color on line) (a) The luminescence, visible and AFM images of XRL beam focusing spots on the surface of a LiF 
crystal obtained after one and three shots with XRL laser intensity 5 mJ/cm ^. Trace was done for AFM image, obtained in 
the case of three shot irradiation of a LiF crystal. 



istry of Education, Science, Sports and Culture, Grant- 
in- Aid for Kiban A No 20244065, Kiban B No. 21360364 
and by the RFBR grant No. 09-02-92482-MNKS-a (AYaF, 
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